Voltage controlled spin injection in a (Ga,Mn) As/(Al,Ga) As Zener diode 
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The spin polarization of the electron current in a p-(Ga,Mn)As-n-(Al,Ga) As-Zener tunnel diode, 
which is embedded in a light-emitting diode, has been studied theoretically. A series of self-consistent 
simulations determines the charge distribution, the band bending, and the current-voltage charac- 
teristics for the entire structure. An empirical tight-binding model, together with the Landauer- 
Biittiker theory of coherent transport has been developed to study the current spin polarization. 
This dual approach allows to explain the experimentally observed high magnitude and strong bias 
dependence of the current spin polarization. 

PACS numbers: 75.50.Pp, 72.25.Hg, 73.40.Gk 



Spin injection is one of the target applications of fer- 
romagnetic semiconductors, which can serve as a natu- 
ral supply of highly spin-polarized carriers. In partic- 
ular, the most intensively studied ferromagnetic semi- 
conductor (Ga,Mn)As can be grown epitaxially on GaAs 
and the carrier-mediated ferromagnetism in this mate- 
rial provides an elegant way to control the ferromag- 
netic properties by tuning the hole concentration^ In 
the design of spintronic devices, the p-type character of 
(Ga,Mn)As introduces a disadvantage due the low hole 
spin lifetimes in GaAs. It has been shown that inter- 
band (Zener) tunneling from valence band electrons of 
(Ga,Mn)As to the conduction band of GaAs is a way 
to circumvent this disadvantaged^ Recently, an injected 
spin polarization up to 80% at 4.6 K has been demon- 
strated in a (Ga,Mn)As based spin-light emitting diode 
(LED) using Zener tunneling^ Moreover, the degree of 
injected spin polarization exhibits a strong dependence 
on the applied bias. The spin polarization reaches its 
maximum just above the electroluminescence threshold 
and decreases dramatically at higher bias. This effect is 
very interesting for spintronic applications, since it pro- 
vides a manner to control the degree of injected spin po- 
larization with the applied voltage. 

In this Rapid Communication, we analyze theoretically 
the transport in the spin-LED as a function of the ap- 
plied bias by means of self-consistent simulations. Fur- 
thermore, we compute the degree of current spin po- 
larization at the (Ga,Mn)As/GaAs-interface by combin- 
ing an empirical tight-binding model with the Landauer- 
Biittiker theory of coherent transport. We show that 
the decrease of the polarization with bias is caused by an 
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FIG. 1: The degree of injected spin polarization measured 
at 4.6 K as a function of the applied bias voltage (a) and the 
current (b) from Ref. 0- 



increased electron tunneling from the valence band of de- 
pleted (Ga,Mn)As and non-magnetic GaAs. This offers 
the opportunity to tune the spin current polarization by 
an external electric field. 

The device considered here has the following struc- 
ture: p + GaAs substrate / 200 nm p-Alo.3Gao.7As (2 x 
10 18 cm- 3 ) / 100 nm p-GaAs (2 x 10 18 cm" 3 ) / 60 nm 
n-ALjGai-zAs (1 x 10 17 cm -3 ) / 30 nm n-Al x Gai_ a As 
(1 x 10 18 cm- 3 ) / 9 nm n-GaAs (9 x 10 18 cm- 3 W 20 nm 
Gao.92Mno.08As, i.e., the spin-LED from Ref . [J. In this 
structure, the Al-concentration in the spin-drift region is 
engineered together with the doping concentration in or- 
der to provide an effective barrier for the holes, such that 
carrier generation due to impact ionization is eliminated 
at low bias. The measured* spin polarization is shown 
as a function of the applied bias voltage and current in 
Fig. n^a) and fl^b). respectively. The spin polarization 
decreases dramatically with increasing bias. 
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FIG. 2: The measured current (solid line) at 4.6 K as a func- 
tion of the applied bias voltage (Ref. 0) and the simulated 
current (dashed line) as a function of the applied bias volt- 
age. The insets show the simulated conduction and valence 
band profiles at (a) 1.8 and (b) 3.0 V bias. The grey-colored 
area represents the (Ga,Mn)As injector. 



We first model the charge transport in the device by 
means of a series of self-consistent simulations of the 
entire LED under bias using Medicij^ a semiconductor 
simulation tool that allows self-consistent calculations 
of semiconductor heterostructures taking into account 
band-to-band tunneling, recombination and impact ion- 
ization. In these simulations the (Ga,Mn)As region is 
treated as a heavily doped GaAs region. The value of 
the doping level is chosen based on transport measure- 
ments on samples grown using similar growth conditions 
on a semi-insulating substrate. In the self-consistent 
simulations, the interband tunneling is taken into ac- 
count using the Kane model to introduce a generation 
ternipi In the simulations, the generation of carriers due 
to interband tunneling is calculated assuming a com- 
pletely filled valence band. Due to the large Fermi en- 
ergy (Ep) in (Ga,Mn)As, the carriers depart from Ep 
below the valence band maximum. This results in a 
small error in the bias voltage where interband tunnel- 
ing takes place, but does not qualitatively alter the re- 
sults of the calculations. Impact ionization is treated in a 
post-processing mode, i.e., the carrier generation is cal- 
culated after the self-consistent calculation of the band 
bending. For comparison, the simulated current density 
at 40 K is shown together with the measured current den- 
sity at 5 K as a function of the applied bias voltage on 
Fig. [3 Due to well-known numerical convergence issues, 
we have been unable to simulate the structure at 5 K. 
There is a good qualitative agreement between the sim- 
ulated and the measured behavior, taking into account 
that the doping levels and aluminum concentrations used 
in the simulations are nominal and can differ from real- 
ity. Between 1.6 V and 2.3-2.5 V the current is dom- 
inated by interband tunneling of electrons through the 
(Ga,Mn)As/GaAs tunnel diode, while above 2.6 V the 
current is dominated by the holes, which can freely flow 
from the substrate to the top contact. Figure|3shows the 
band diagram (solid lines) of the (Ga,Mn)As-(Al,Ga)As 
diode together with the Al concentration and the car- 
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FIG. 3: The Al concentration (solid line) (a), the band di- 
agram (solid lines) and the carrier generation (filled) due to 
interband tunneling (dotted line) as a function of the depth 
near the (Ga,Mn)As/GaAs tunnel diode at 1.8 V (b) and 3 V 
(c). The grey-coloured area represents the (Ga,Mn)As in- 
jector. The dotted arrows indicate the tunneling of valence 
electrons from (Ga,Mn)As to GaAs. 



rier generation due to interband tunneling (filled black) 
as a function of the depth near the (Ga,Mn)As/GaAs- 
interface at 1.8 V and at 3 V. There are two peaks present 
in the plotted carrier generation, one in the (Ga,Mn)As- 
area, which shows the generated holes in the (Ga,Mn)As 
valence band, and one in the GaAs area that shows the 
generated electrons in the GaAs conduction band. From 
this picture we can deduce that interband tunneling takes 
place from the valence band of (Ga,Mn)As to the conduc- 
tion band of GaAs. When GaAs changes into (Al,Ga)As, 
the increasing Al concentration causes a widening of the 
band gap and an abrupt increase of the tunneling dis- 
tance. This results in an exponentially smaller tunnel- 
ing probability and hence the number of electrons that 
is generated in the (Al,Ga)As region is negligible. This 
means that only the part of the (Ga,Mn) As valence band 
that overlaps with the GaAs conduction band can par- 
ticipate in the tunneling process. If we compare the low 
and the high bias case we see that the area from where 
electrons tunnel (the left peak in the carrier generation) 
is much bigger at high bias than at low bias. When the 
voltage drop over the tunnel diode increases, the part of 
the valence band of (Ga,Mn)As that aligns with the con- 
duction band of GaAs increases and hence the tunneling 
originates from a wider region. In Fig. 0] we show the 
hole concentration and the (normalized) carrier genera- 
tion due to interband tunneling as a function of the depth 
near the (Ga,Mn)As-GaAs interface for different values of 
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FIG. 4: The hole concentration (left) and the carrier genera- 
tion due to interband tunneling (right) at different values of 
the applied bias voltage near the (Ga,Mn)As/GaAs interface. 
The grey-coloured area represents the (Ga,Mn)As injector. 



the applied bias voltage. The simulations show that the 
hole concentration decreases in the bottom two nanome- 
ter of the (Ga,Mn)As-layer due to the band bending in 
the (Ga,Mn)As-GaAs p-n-diode. Increasing the voltage 
has a negligible influence on the hole depletion near the 
interface. However, we do see a big change in the prop- 
erties of the tunneling region. At low bias, the tunneling 
electrons mainly depart from a region where the hole con- 
centration is smaller than, but close to the bulk value, 
while at higher bias this region shifts and widens such 
that electrons can tunnel from a region where the hole 
concentration is much lower. At a somewhat higher bias 
(> 2.2V) we see that also electrons from the valence band 
of the non-magnetic GaAs layer participate in the tun- 
neling process. In the experiment this bias will probably 
be lower due to the above described small error in the 
interband tunneling model. 

Apart from Zener tunneling, also impact ionization 
can contribute to the generation of electrons in the 
GaAs conduction band. This process generates un- 
polarized electron-hole pairs and hence dilutes the in- 
jected spins and diminishes the measured spin polariza- 
tion. Impact ionization in this case is caused by holes 
that flow from the substrate to the top contact and are 
heavily accelerated by the strong electric field near the 
(Ga,Mn)As-GaAs interface. However, the simulations 
indicate that the carrier generation due to impact ion- 
ization only starts dominating at the (second) hump in 
the IV-characteristics, where the valence band reaches a 
flatband-situation (Fig. |3 inset (b)) and the holes can 
flow freely from the substrate to the top contact. Below 
this "hump" impact ionization can be neglected. 

In parallel to self-consistent calculations that allow a 
detailed understanding of the charge distribution, band 
bending, electron and hole currents and their bias de- 
pendence, in order to study the current spin polariza- 
tion, we have developed a model of interband tunnel- 
ing, which combines an empirical tight-binding approach 
with the Landauer-Biittiker formalism. 8 To describe the 
band structure of GaAs we use the sp 3 d 5 s* tight-binding 
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FIG. 5: Scheme of the structure used in the tight-binding 
calculations, at zero and at reverse bias without (a) and with 
the GaAs spacer (b). 



parametrization, with the spin-orbit coupling included, 
proposed by Jancu et al. 9 This model reproduces cor- 
rectly the effective masses and the band structure of 
GaAs in the whole Brillouin zone, in agreement with the 
results obtained by empirically corrected pseudopotential 
method. The presence of Mn ions in GaMnAs is taken 
into account by including the sp-d exchange interactions 
within the virtual crystal and mean-field approximations, 
with the values of the exchange constants determined by 
the observed spin splittings of the conduction and va- 
lence bands. Importantly, in contrast to the standard 
k -p method, our model takes automatically into account 
the Rashba and Dresselhaus terms, and is therefore par- 
ticularly well suited to describe interface phenomena such 
as tunneling. 

Because of computational constraints in approaches 
involving transfer matrix formalism, 8 carrier transport 
along the whole device cannot be simulated. Therefore, 
we consider first the simplest p-Gai-^Mn^As/n-GaAs 
tunneling structure shown in Fig. ISJa) . Guided by previ- 
ous theories of tunneling magnetoresistance (TMR), we 
expect that while such model overestimates necessarily 
the tunneling current, it can provide quantitative infor- 
mation on current spin polarization that is determined 
by inter-band tunneling matrix elements and degree of 
spin polarization in the ferromagnetic electrode. 

In our computation, we assume that the magnetization 
vector is by 27° out of plane, as implied by experimen- 
tal conditions^ We evaluate the spin current polarization 
Pj in respect to this direction. Furthermore, we estimate 
that for a T c of 120 K and x — 0.08, the expected hole 
concentration is of the order of p = 3.5 x 10 20 cm -3 , as 

We also 
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indicated by the experimental results in Ref. 10 
assume that the electron concentration is n = 10'" cm 
The dependence of Pj on the hole concentration p is de- 
picted in Fig. |SJ Previous calculations have shown that 
the spin polarization of the hole liquid decreases with p 
in Gai-^Mn^Asii In this case, Pj, the spin polarization 
of the tunneling electrons, shows a similar behavior as 
function of p. This is naturally caused by the fact that 
in the relevant range of p and x, Gai-^Mn^As is not 
half-metallic, so that the Fermi energy is greater than 
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FIG. 6: The calculated bias dependence of the current polar- 
ization for different hole concentrations in the Gai-^Mn^ As- 
layer for structure (a) of Fig. 



the spin-splitting. Surprisingly, however, a weak but op- 
posite behavior is seen at low Vz- We assign this result to 
the fact that the current polarization is determined not 
only by the electron spin polarization at the Fermi energy 
but also by the selection rules for transition probabilities. 
Actually, we know that the mixing of the spin wave func- 
tions, which reduces the tunneling matrix elements, is 
only important if the Fermi energy is much smaller than 
the spin-orbit splitting of the valence band. Accordingly, 
the reduction of Pj by the spin-orbit coupling decreases 
gradually with the hole density. Nevertheless, this effect 
appears to be too weak to explain much smaller Pj in the 
non-annealed sample from Ref. 0. 

Next, to simulate more realistically the device, we con- 
sider the structure consisting of Gai-^Mn^As/GaAs/n- 
GaAs, where the GaAs spacer has to be kept thinner 
than that of the n-GaAs depleted layer implied by our 
self-consistent calculations. The schematic of such struc- 
ture is shown in Fig.[S{b). Figurc[7|shows Pj as a function 
of the bias Vz of the Zener tunneling diode for various 
thicknesses d of the spacer layer. We see that at low bias, 
Pj depends weakly on d and is of the order of 0.7, in per- 
fect agreement with the experimental results. Interest- 
ingly, the drop of Pj with Vz , is greater for larger d. The 
strong dependence of Pj on Vz is again in agreement with 
the experimental findings, though a direct comparison is 



hampered by the fact that the exact relation between Vz 
and the total bias V applied to the device is unknown. 
The dependencies of Pj on Vz and d can be easily ex- 
plained with the help of Fig.[5Jb), which shows that for 
Vz > the holes tunnel partly from the non-magnetic 
GaAs. 

In conclusion, we have analyzed the spin polariza- 
tion of the injected current in a (Ga,Mn)As-(Al,Ga)As 
spin-LED by performing self-consistent simulations of the 
band bending and using tight-binding model together 
with the Landauer-Biittiker formalism for the calcula- 
tion of the tunneling current. Our studies explain quan- 
titatively the large spin polarization of the injected cur- 
rent and its strong dependence on the applied voltage, 
1.0 



o 0.8 

"5 



o 

D. 

E 0.4 

CD 



o 



0.2 



0.0 



T 



T 



1 I 
monolayers 
8 monolayers 
16 monolayers 



0.0 



0.1 



0.2 



0.3 



0.4 



0.5 



V Z (V) 



FIG. 7: The calculated bias dependence of the current polar- 
ization for different thicknesses of the spacer layer, for struc- 
ture (b) of Fig. U 



as has been recently observed in experiments. We as- 
cribe the observed decrease of the spin polarization with 
the applied bias to the enhanced tunneling from depleted 
(Ga,Mn)As and non-magnetic GaAs regions. This can 
provide a path towards voltage controlled magnetic be- 
havior. By choosing the right voltage drop over the tun- 
nel diode, one can switch between injection from a mag- 
netized region to injection from a non-magnetized region. 
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